Des ondes continues aux matrices discrètes
La vibration d'une corde ou d'une membrane est régie par l'équation des ondes :
$$\rho(x) \frac{\partial^2 v}{\partial t^2}(x, t) = \frac{\partial}{\partial x} \left[ p(x) \frac{\partial v}{\partial x}(x, t) \right]$$
Pour trouver la solution $v(x, t) = \sum_{k=0}^{\infty} c_k u_k(x) \cos \sqrt{\lambda_k}(t - t_0)$, nous devons résoudre le système de Sturm-Liouville:
$$\frac{d}{dx} \left[ p(x) \frac{du_k}{dx}(x) \right] + \lambda_k \rho(x) u_k(x) = 0$$
La discrétisation de l'opérateur conduit à des équations matricielles telles que $Aw = -0.04 \frac{\rho}{p} \lambda w$. Pour une matrice tridiagonale $4 \times 4$, $p(\lambda)$ est gérable. Cependant, lorsque le maillage se raffine ($n$ augmente), nous rencontrons deux obstacles :
- Limite de Abel-Ruffini : Aucune solution algébrique n'existe pour les racines des polynômes dont le degré est $n \ge 5$.
- Sensibilité aux arrondis : Dans les systèmes de grande dimension, un changement dans la dixième décimale d'un élément peut déplacer les valeurs propres d'un ordre de grandeur (phénomène du polynôme de Wilkinson).
Nécessité numérique et bibliothèques professionnelles
Les bibliothèques numériques professionnelles (IMSL, NAG) évitent les polynômes caractéristiques bruts. Elles utilisent plutôt des algorithmes itératifs pour l'approximation :
- Bibliothèque IMSL : Utilise la régression linéaire au moindre carrés, les splines cubiques et les transformations de Fourier rapides.
- Bibliothèque NAG : Emploie l'approximation polynomiale au moindre carrés et les ajustements selon les normes $l_1/l_{\infty}$.
Lors de l'approximation des valeurs propres du système $\lambda_i = 1 + 4\alpha\left(\sin \frac{\pi i}{2m}\right)^2$, nous nous appuyons sur les moindres carrés discrets et la découverte itérative plutôt que sur la recherche de racines.